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Abstract 



A pre-synaptic neuron releases diffusing neurotransmitters such as gluta- 
mate that activate post-synaptic receptors. The amplitude of the post- 
synaptic current, mostly mediated by glutamatergic (AMPARs) receptors, 
is a fundamental signal that may generate an action potential. However, al- 
though various simulation results 
control the post-synaptic current, it is still unclear how this current de- 
pends analytically on factors such as the synaptic cleft geometry, the dis- 
tribution, the number and the multi-conductance state of receptors, the 
geometry of post-synaptic density (PSD) and the neurotransmitter release 
location. To estimate the synaptic current maximal amplitude, we present a 
semi-analytical model of glutamate diffusing in the synaptic cleft. We mod- 
eled receptors as multi-conductance channels and we find that PSD mor- 
phological changes can significantly modulate the synaptic current, which 
is maximally reliable (the coefficient of variation is minimal) for an optimal 
size of the PSD, that depends on the vesicular release active zone. The exis- 
tence of an optimal PSD size is related to nonlinear phenomena such as the 
multi-binding cooperativity of the neurotransmitter to the receptors. We 
conclude that changes in the PSD geometry can sustain a form of synaptic 
plasticity, independent of a change in the number of receptors. 
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Introduction 

Synapses are local active micro-contacts underlying direct neuronal com- 
munication. Depending on the brain area and the neuron types, synapses 
can vary in size and molecular composition. These inter-synaptic varia- 
tions are mediated by hundreds of different molecules and proteins, partic- 
ipating in the assembly of the stable but plastic synaptic structure 
Neurotransmitters such as glutamate molecules, after being released from 
vesicles, diffuse in the synaptic cleft, between the pre and post synaptic 
terminals (FigfT]). The post-synaptic terminal of excitatory synapses con- 
tains ionotropic receptors such as AMPA and NMDA receptors and they 
may open upon binding with neurotransmitters. AMPARs are tetrameric 
assemblies composed of four different subunits, which can bind to a gluta- 
mate molecule (0), but it has been reported that two agonist molecules at 
least are required to open a sing le AMPA channel fl). The amplitude of 
ionic current is thus proportional to the number of open receptors and their 
conductances. The postsynaptic current measures the efficiency of synaptic 
transmission and reports in a complex manner the frequency and location 
of released vesicles p). 

It is intriguing that although the number of neurotransmitters released is 
of the order of thousands, the number receptors is at most one hundred. This 
difference may serve to compensate the small receptor patches that should 
be found by the neurotransmitters 113) • To study synaptic transmission, a 
fundamental step was the analysis of channels such as AMPARs, expressed in 
oocytes. By recording the current using patch clamp and excised patch, the 
conductance state properties, related to the open, closed and desensitized 
states have been extracted by Markov chain models (@-@). However, to 
reduce the complexity of AMPARs dynamics, Markov chains were kept as 
minimal as possible, based on one or two possible binding sites (0-0) with 
a single conductivity level. Only recently (|10l ). to study the high variability 
(~ 5 — 100 pA) of the synaptic current I s , a four state channel model was 
used to interpret the high coefficient of variation, which was shown to be 
due to the spatiotemporal correlations of two released vesicles. When the 
receptor properties have been sufficiently well characterized, a second step 
consisted in integrating these properties within the synaptic organization to 
reconstruct the synaptic function. This step became possible by the use of 
modeling and numerical simulations 
role of synaptic geometry on the number of open receptors. More recently 
using glutamate diffusion and electrical resistance properties of the synaptic 
cleft, it was anticipated that the synaptic current could be maximal for an 
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optimal cleft height (|2C 

We find here by deriving a novel semi-analytical approach that the synap- 
tic current depends on the vesicular release location, the number and the 
biophysical properties of the receptors, the PSD size and location, and the 
geometrical characteristics of the synaptic cleft. To obtain these estimates, 
we model glutamate as diffusing molecules and approximate the cleft geom- 
etry as a narrow cylinder (j2ll ). But instead of using the classical Markov 
description (0-H), our analysis relies on some direct analysis of AMPA con- 
ductances states @) and we account for the four glutamate binding sites 
per receptor. Our analysis reveals that given the pre-synaptic active zone 
size, where vesicles are released, the coefficient of variation (CV which is the 
standard deviation divided by the mean) of the synaptic current is minimal 
for a specific PSD size and all other sizes, that may be induced by plastic 
changes (|23|), lead to an increasing CV. However, for a centered active zone, 
we show that for a fixed density of receptors, the synaptic current is always 
a decreasing function (to zero) of the PSD radius. We further show that a 
maximal and reliable current cannot be achieved simultaneously for the same 
distribution of synaptic parameters. Finally, we will propose that a synapse 
can increase its reliability by restricting the active zone (AZ) radius, which 
is a unique and nonlinear function of the PSD radius. This result should 
be true for generalized geometry and not only geodesic disks. Finally, re- 
modeling the PSD, which affects the synaptic current, can occur in parallel 
with the classical synaptic modulation, induced by the direct addition or 
removal of synaptic receptors. PSD remodeling can be mediated by geomet- 
rical or internal scaffolding reorganization that can be transient but 
much faster than changing permanently the synaptic receptors. These fast 
changes can thus affect the detection threshold of the post-synaptic neuron 



24 . l25l ) and be a source of synaptic plasticity without changing the 



number of AMPARs. 



Method: Theoretical model 
Diffusion in the synaptic cleft. 

The synaptic current I s is mediated by open AMPARs, which can bind 
from one to four glutamate molecules. A single AMPAR has several con- 
ductance states which correspond to the combination of the four distinct 
conductances, associated to the different GluR subunits (5 . l26h. However, 
to reduce the complexity of the analysis, most studies da-la, llll, uA, Il5l ) have 
modeled AMPARs dynamics by one or two bound glutamate molecules. The 
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conductances designated by 71,72,73,74 associated with 1, 2, 3 or 4 bound 



glutamate molecules. It has been reported (|27l . |28| ) that at least two glu- 
tamate molecules are needed to open an AMPAR, so we set 71 = 0. The 
synaptic current I s (t) depends on the number of open channels (./V2, A^, N4) 
bound respectively by 2,3 or 4 glutamate molecules. In that case, 

h(t) = (j 2 N 2 (t) + j 3 N 3 (t) + l4 N A {t))AV, (1) 

where AV is the difference of potential between the intra and the extracellu- 
lar medium. We selected from the experimental results (0), obtained for dif- 
ferent glutamate concentration the values 72 = 4p5, 73 = WpS, 74 = 13p5 
(for AV m = -lOOmV). 

Our goal here is to estimate the mean and the variance of open receptors 
and quantify the peak amplitude of the current I s . After a vesicle fuses 
with the pre-synaptic membrane at position xq, N g = 3000 glutamates are 
released (FigfT]). Glutamate molecules diffuse in the cleft and are reflected 
on the synaptic membrane, while they are absorbed at the lateral boundary 
of the synaptic cleft. AMPARs are uniformly distributed over the PSD. We 
consider that only a certain fraction of glutamates hitting an AMPAR leads 
to receptor activation, due to a chemical energy barrier. The synaptic cleft 
is modeled as a cylinder £1 (FigfT]), while there are N a AMPARs located 
on the PSD. When a glutamate molecule hits an AMPAR, it can either be 
reflected or it will activated the receptor. We model this behavior using 
a homogenized radiative boundary condition over the entire PSD (OQpsd) 
A glutamate that hits the neuron membrane is reflected except at 



the lateral cleft boundary (dtliat), where it will not contribute to activate 
an AMPAR and this is modeled by an absorbing boundary condition. 

Estimating the number of open receptors. 

The probability to find a glutamate molecule at position x at time t, when 
it started at position xq is given by the density function p(x,t\xo) that 
satisfies the equation: 



dp(x,t\x ) 



dt 

p(x, 0) = 5(x — xq 



DAp(x,t\x ), x e t > (2) 



dp{x,t\x ) 



-D 



dv 

dp(x,t\xo) 



dv 



0, p(aj,t|aJo)lan iat = 

-Kp(x,t\x ), 
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where D is the free glutamate diffusion constant. The partial absorption 
constant k accounts for the fraction of AMPARs inside the PSD and the 
activation barrier of a glutamate to a GluR binding site. Using a homog- 



enization procedure (|3CH 32|. |4Q| ). we shall derive a new expression in the 

(3) 



context of the synaptic cleft for the partial reflection parameter k 

D 1 



ybL> N a a ^ Ka2TTa 2 N a 

where f(a) = 1 — a, a = N a a 2 / R 2 PSD , a is the radius of the binding site 
and Rpsd is the radius of the PSD and K a measures the partial binding of a 
glutamate molecule to a single receptor. We have summarize in the appendix 
all these steps. To determine the constant «„, we used the fitted results 
obtained from the Markov analysis given in (|7j). In a first approximation 
neglecting n a , the partial reflecting constant k is directly proportional to 
the binding rate of glutamate molecules and in the appendix, we obtain the 
numerical approximation K a ~ 1.06. 

To determine the value of k (see appendix 1 1 . 2|) we further need to esti- 
mate the effective radius a of a single receptor. For that purpose, we run 
some simulations for a typical synapse of radius 500nm with a PSD radius 
of 300nm, a height of 30nm and an AZ radius of 150nm. Using the criteria 
that the synaptic current saturates for four released vesicles (~ 12, 000 glu- 
tamate molecules), we obtain that radius a ~ 1.8nm, as shown in figure ([3]). 
Interestingly, the radius a accounts not only for the geometrical properties 
of the AMPAR binding site, but also for the underlying electrostatic inter- 
actions. This value a should be compared to the recent crystal structure 
dimensions (the AMPAR has a transversal size of 9nm and the total length 
is around 18nm) of the ligand-binding domain reported to be less than 4nm 

B). 

The probability p(xq) that a glutamate molecule released at position xo, 
binds a receptor is given by the total flux: 

JO JdUpsD 0n JdQpsD 0n 
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where u(x\xq) = J °° p(x, t\xo)dt satisfies 



DAu(x\xq) 



-5(x — xq) for x G 



du(x\xo) 



ch 



0. 



u(ce|;eo) 



an r 



D 



du(x\xo) 



dv 



-KU{X\Xq), 



dflpsi 



In the appendix, we derive an analytical expressions of u and p. For a vesicle 
releasing its contain at the center of the AZ, the probability to bind one of 
the AMPARs is 



P 



JPSD 



JpSD + Jhat 



where 



and 



Jpsd = 2-7TK / -^-K (ar)rdr + 2ttk f AI ( 
Jo D Jo 



ar)rdr 



Jl 



at 



-2-nDhC. 



where A and C are two constants that depend on L, k, D, and h, given in 
the appendix. 



The mean and variance of the synaptic current I s . 

To compute the mean and the variance of the maximal amplitude of synaptic 
current I s , we shall account for two possible sources of fluctuations: one is 
due to the number of bound glutamate molecules and the second to the con- 
figuration of bound AMPARs. By configuration, we mean the distribution 
of AMPARs bound to 2, 3 and 4 glutamate molecules. 

To estimate the first source, we use the probability distribution Prk{xo) 
to have k glutamate molecules bound, when a vesicle is released at position 
Xq. This probability follows a binomial distribution 

Pr k (x ) = C k Ng p(x ) k (l-p(x )) N s- k (4) 



p{x ) = k u(x\x )dS{x), (5) 

JdQpsD 
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where N g is the total number of released glutamate molecules. Because 
AMPARs can bind from zero to four glutamate molecules, the probability 
of a given configuration n = (714, 7^3, 712, t&i) to have 71 1 AMPARs bound to 
one glutamate, 712 AMPARs bound to two glutamates and so on, when there 
are N a AMPA receptors, for k < min(4N a , N g ) bound glutamate molecules, 
is given by 

N 1 1 

Pr-fnlAfl- = — — (q) 

111 n 4 Wn 2 \n 1 \(N a -(n 4 + n 3 + n 2 +n 1 ))\F(k,N a y K) 

where this probability is computed by choosing n 4 AMPARs out of N a , n 3 
out of N a — n 4 and so on. F(k, N a ) is the number of possibilities to decompose 
the integer k on the integer 4,3,2, 1,0, when there are at most N a terms: 
that is k = 4n4+3ri3+2rt2+7ii. In practice, we compute F(k, N a ) numerically 
as the k + l's coefficient of the expression (1 + x + x 2 + x 3 + x 4 ) Na . The 
present analysis can be used to obtain any statistical moments associated 
to the current and in particular, the mean and variance 

(J s (aj )> = n-7PT{n\k}Pr k (x )= (7) 

k=l n€Sk 
4N a 

AV^2 S2 n ■ jPr{n\k}Pr k (x ) + 

k=l nS5fe 



/ 47V a \ 

AVN al4 \l-J2 p rk(xo)\ 

V k=0 J 



{I 2 s (x )) = AV 2 ^J2^-^ 2pr i H \ k } Pr k( x o)-(Is(^)f = 

k=l n&S k 

N a 

J^Y,^- ^ Ay ) 2 Mn\k}Pr k (x ) + 

k=l n£Sk 

( 4Na \ 

(N al4 AV) 2 1 - Y,Pn(x ) - (Is(x )) 2 

V k=0 J 

where AV is the voltage drop, 7 = (71, 72, 73, 74) is the conductances vector, 
Sk represents the set of possible configurations of n = (711,712,713,714) such 
that 4n 4 + 3rt3 + 2n2 + n\ = k. The formulas for the mean and variance are 
made of two terms: the first is the sum over all the sites that are partially 



Multi-conductance AMPA receptor model 



7 



bound by glutamate molecules and the probability for such an event is the 
product of the probability Pr^ that k glutamates are bound (k < 4A r a ) and 
the probability Pr{?i|/c} of a given binding configuration k = An^ + 3n^ + 
2n2 + n\. The second term accounts for all bound AMPARs (4iV a ) and this 
happens with the complementary probability to the first case. 

Numerical simulation results for the synaptic cur- 
rent 

To determine the maximum amplitude of the synaptic current, we use our 
semi-analytical method developed above and simulate a single vesicle (N g = 
3000) released. We find (FigJ2R) that AMPARs are mostly bound to two 
glutamate molecules, while for already two vesicles, the dominant contribu- 
tion comes from receptors bound to four. To further investigate the influ- 
ence of the PSD size on the synaptic current, We plotted the current I s as 
a function of the release site position (Figj2j3): when release occurs outside 
the region above the PSD, the current drops drastically due to a fast decay 
of the binding probability. 

Since our analysis allows us to determine the relative influence of the 
AZ and the PSD, we vary their respective sizes and we estimate the conse- 
quences on the synaptic current. We first plotted in figureEA, the mean and 
variance of I s as a function of the PSD radius when one vesicle is released 
at the center. We find that for a given density of receptors, the current is a 
decreasing function of the PSD radius. In figure [4)3-C, we plotted the num- 
ber of AMPARs bound to two, three and four glutamate molecules when we 
fixes the AZ radius to 50nm (blue) and 150nm (red). For a small AZ radius, 
the AMPARs are in majority bound to four glutamate molecules and thus 
the synaptic current amplitude is much higher compared to the case of a 
large AZ radius. In that latter case, the current is primarily generated by 
receptors bound to two glutamate molecules. We can now use our refined 
analysis to re-interpret the large current differences observed in figure HA 
with the radii 50nm and 150nm. This difference is due to the nonlinear prop- 
erties of having different conductivities, generated by the amount of bound 
glutamate molecules. But in all cases, the synaptic current is a decreasing 
function of AZ radius (Fig. HP). 

To asses the reliability of the synaptic response, we use the coefficient of 
variation (CV=standard deviation over the mean) of the synaptic I s . For a 
fixed AZ radius, we now show that the CV has a minimum as function of the 
PSD radius: Indeed, using our simulation results HJC, we find for example 



Multi-conductance AMPA receptor model 



8 



with an AZ of radius 100 nm, that the CV reaches its minimum for a PSD 
radius of 120 nm. In order to confirm this result, we performed Brownian 
simulations associated with our semi-analytical model. This Brownian sim- 
ulation is described as followed: glutamate molecules are released from the 
AZ of radius 250nm and diffuse in a synapse of radius 500nm. We run sim- 
ulations (figHJK) for various PSD sizes and computed the associated CV (for 
an average of 100 samples). In the simulations, a receptor is modeled as a 
partial absorber and when it binds to four glutamate molecules, it becomes 
totally reflecting. The current is computed using equation [TJ Finally, any 
glutamate molecule reaching the lateral boundary of the synapse is perma- 
nently absorbed. Using these rules, We find (figj6f3) that our initial analysis 
is confirmed by the Brownian simulations . 

To further investigate the relation between the radius of the PSD and 
the AZ, we plotted (figure UP) the optimal PSD radius as a function of 
the AZ size. We find that the optimal PSD size increases with the size of 
the AZ, but the relationship is a nontrivial nonlinear. To see whether the 
size of a synapse affects the CV curve, we fixed the AZ radius at 50 nm 
and plotted the synaptic current as a function of the PSD radius for four 
different synapses of sizes 200,300,400 and 500 nm (Figj5|). Interestingly, as 
shown in Figj5j>, the value of the PSD for which the optimal CV is achieved, 
does not depend much on the synaptic radius. We observe that the CV is a 
decreasing function of the synaptic radius and thus large synapses are more 
reliable than small ones. Finally, to assess the role of a possible fluctuation 
in the number of glutamate molecules released from a vesicle, we show in 
figure Figj5p the effect of three different distributions for 2000, 3000 and 
4000. This shows very little change in the CV minimum phenomena. In 
Figj5p, we compare release where the mean number is 3000 molecules with 
a variance of 500 Gaussian distributed (bold line) with a release of a fixed 
number 3000 (dashed line). We observe a small deviation, showing that in 
this range of fluctuation, the number of released glutamate molecules does 
not affect much the CV. 



Discussion 



We have analyzed here the synaptic current I s starting from the intrinsic 
biophysical properties of the AMPARs (0). We included the diffusion of 
glutamate molecules in the narrow synaptic cleft and estimated the fraction 
that binds AMPARs. Contrary to previous works (0, El , 14, 15), our analysis 
does not use the description of channels involving a time dependent multi- 
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states Markov chain, instead we use a multi-conductance approach mixed 
with a time independent receptor description. In addition, we neglected any 
possible interactions between bound receptor subunits that would affect the 
probability for a free subunit to bind a glutamate molecule. This approach 
allows us to obtain semi-analytical results about the synaptic current and 
to account for the nonlinear effect due to the multi-conductance states of 
bound AMPARs. We also studied the role of the cleft geometry and explore 
the role of several parameters. 

For example, we have found here that modulating the PSD size can 
affect the synaptic current amplitude: the current depends on the relative 
size of the PSD to the AZ (Fig. 0]A.,D). As a consequence, because the PSD 
could be reorganized, while the total number of receptors remains constant, 
we propose that this direct change in the PSD geometry can modulate the 
synaptic current. These functional consequences were anticipated in recent 
experiments using GFP-tagged PSD-95 (1231 ). As suggested there, the PSD 
size is in constant remodeling. Thus, combined with our analysis, we propose 
that synaptic changes are a source of synaptic current modulation. It is not 
clear what are the reasons for these changes, but they could be induced by 
long term potentiation (LTP) or depression protocols (|23l). Changing PSD 
size, while keeping AZ fixed can be seen as a form of plasticity induced by 
structural remodeling without any change in the number of receptors. These 
changes can, for example, be induced by actin dynamics, which is correlated 
with spine shape changes (|34l . l35l ) . 

We have not accounted here neither for glial cells nor neuronal trans- 
porters, that only weakly affect direct synaptic transmission and the num- 
ber of open AMPARs. The effect is of the order 10% - 15% ((0),(0) 
Figure 9, (0), ( 17). (fl9l )). However, in some pathological conditions, related 
to a glial reorganization, glial transporters can directly modulate synaptic 
transmission (J9j), which changes drastically the present results. 

Another aspect of our modeling approach relies on the description of AM- 
PARs. Indeed, we computed the synaptic current from conductances orig- 
inating from patch-clamp experiments of isolated AMPARs (0, 129] ). where 
the relation between the number of bound glutamates and the associated 
conductances was obtained for different fixed glutamate concentrations. Al- 
though more than four conductance levels have been reported, it is still 
unclear how to relate them to the number of bound glutamate molecules. 
In that context, it would be interesting to design a specific experiment to 
measure simultaneously on single AMPAR the number of bound glutamate 
and the associated current. In addition, having four glutamate bound to 
a single receptor lead to an amplitude of 13 p A, which has to be compared 
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with two bound AMPARs to two glutamates leading to 2 • 4 = 8pA). This 
difference suggests that binding four glutamates to a single AMPAR has a 
nonlinear effect and is different from having two receptors bound by two glu- 
tamates, as reported in figures H]A[3J3 and SC. To show that the minimum 
of the CV comes from this nonlinearity, we shall now consider a reduced 
model for single receptor with three states, immersed in an ensemble of N 
binding molecules. Depending on the number of bounds, the receptor can 
be in one of the three states, where the current I can switch between the 
three values [Ii, I2, 13}. Each molecule binds with a probability q 6 [0, 1] and 
the probability of / is given by 



p(I = h)={ 



^\ k {l-q) N ~ k it = 1,2 

) (8) 



[ l-p(I = h)-p(I = h)-(l-q) N k = 3. 

where the probability for one and two is given by the Binomial law while to 
compute the probability for the current ^3, we use that it is the difference 
of the first twos minus the probability that no molecules are bound. We 
computed analytically the CV of the current / as a function of the probabil- 
ity q and we observed (Fig|6j3) that only for some range of the parameters 
such as I\, I2 <C ^3, the CV presented a local minimum. This results show 
this simple model captures the features for a minimum. In addition, in that 
model changing the parameter q is equivalent to vary AZ or the PSD radius. 

To conclude, it is still unclear what defines the detection threshold of 
a post-synaptic neuron. Indeed, this current is mediated by the number 
of open AMPARs and the associated conductances. When an AMPAR is 
maintained at the PSD by scaffolding molecules such as PSD-95, located 
just underneath the location of vesiclar released (a signal that can be medi- 
ated by N-cadherin molecules ({22! )). the probability of glutamate binding is 
maximal and thus the AMPARs will report more accurately this vesicular 
event. We predict that this effect will be increased when the scaffolding 
molecules PSD-95 will be over-expressed. Indeed the over-expression will 
results in an increasing of the number of anchored receptors and their rel- 
ative location in comparison with the pre-synaptic terminal. In contrast, 
in PSD-95 knockdown (by shRNA (j25l )). the detection threshold was found 
to be lower, due to a decrease in the number of AMPARs and a disperse 
AMPAR configuration. 

A synapse is an unreliable device but the variability is reduced when 
the AZ and the PSD are apposed with a precise relationship between their 
radius (Fig. EJ 1 ). However, if the vesicular release is spread over the AZ, 
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receptors over the PSD can detect a vesicular event, but it might be very 
small. The regulation of the PSD size should thus be a fundamental param- 
eter comparable to increasing the AMPAR number, which is the molecular 
basis for robust synaptic plasticity. Finally, from a biophysical point of 
view, the partial reflection constant that have introduced for the simulation 
of the AMPARs reflects both the energy activation barrier of glutamate- 
AMPAR interaction and the probability to find a receptor inside the PSD. 
The concept of alignment of vesicular release domain with the PSD was al- 
ready suggested in (0) without quantitative analysis. Thus changing the 
PSD shape can modulate the synaptic current and thus can be considered 
as a source of plasticity. It would be interesting to analyze the molecular 
mechanisms responsible for such changes (]23f ). 



1 Appendix 

Summary of our methodology approach 

We shall now summarize our methodology and the organization of this ap- 
pendix. Our new method consists in combining analytical computations and 
numerical simulations to estimate the probability that a combination of glu- 
tamate molecules activate an AMPAR, with direct experimental measured 
conductances (jl). The methodology to estimate the distribution of bound 
AMPARs is first to compute the probability that glutamate molecule binds 
a receptor. For this we solve in appendix II .1\ the probability equation and 
calculate the flux of receptors to a given region (the PSD) and the total 
flux through all the synaptic cleft as a function of variable x, position where 
glutamate molecules are released. The ratio of the fluxes is the probabil- 
ity. Second, because the height of the synaptic cleft is small, the solution 
u of equation [9] is close to its average over the z-direction and we use this 
property to perform our computations. Third, we consider that the AM- 
PARs are uniformly distributed over the PSD and we use a homogenization 
procedure to replace by a partial reflecting boundary condition summarized 
in the constant k, a complicated boundary condition, where sometimes a 
glutamate molecule would bind to an AMPAR receptor with a given prob- 
ability of activation and sometimes it would be reflected when hitting parts 
of the PSD containing no receptors. Fourth, the rational expression for k is 
derived in appendix II .21 and the associated mathematical derivation is given 
in appendix 11.51 To check the validity of our computation, we compare our 
results with Brownian simulation, described in appendix 11.31 Finally, we 
use published data to obtain an approximate value for k and relate it to the 
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activation and the effective binding size of a single AMPAR (appendix ll.4p . 

1.1 Analytical expression for the binding probability of a 

glutamate molecule to the PSD, using an averaging method. 

We present here an averaging method to obtain an explicit expression for the 
probability p(xq) that a glutamate binds one of the AMPA receptors before 
it escapes. To compute the total number of bound glutamate molecules to 
the AMPARs, we solve the steady state diffusion equation in the cylindri- 
cal synaptic cleft geometry, where glutamate is released at a^o- Using an 
averaging method, we analyze equation 



DAu(r, z) = -5(r)5{z - z ) for {(r, z)\r £ [0, R), z G (0, h)} (9) 



du(r, z) 



Oi 



0, u(r,z)l 



D 



du(r, z) 



dv 



-Ku(r, z), 



r<L,z=0 

In cylindrical coordinates, the average u(r) = ^ f£ u(r, z)dz satisfies 

1 



1 Id 

D ( u"(r) + -u'(r) - j i Q^u(r,0)x[o,L]{r) 



rh 



6(r) 



(10) 



Using the boundary conditions on dflpsD for r < L, z = 0, we express -§^u(r, 0))x[o,l] ( r ) 
in terms of u: 

u(r, z) « u(r, 0) + du ^°K + {h 2 ). 
oz 

Integrating the Taylor expansion with respect to z, and using that for 
u (r,0) = -g^andr<L, 



,*Mf5 + *|, 



dz \ k 2 
by substituting eq. ([TT|) in (fTTIj) . we get 



fill 



u"{r) + i«'(r) 
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The solution is given by 
u(r) 

where 



2^D K o( ar ) + AI o(ar) < r < L 



r 

C log — L <r < R 

R 



2k 



a 



h(2D + Kh)' 



To determine the parameters A and C, we use the continuity of u and its 
derivative at r = 0. We obtain the linear system to invert: 

^K ( a L) + AIo(aL) = Clog^ 

-—aK^aL) + Aah(aL) = C-. 

The solution to these equation is given by 

A \ ( aI (aL) - log £ V V -±K (aL) 
C ) V ah(aL) x / V 7) aR l( aL ) 

L' 1 K Q (aL) + log ^Kt(aL) 



Da (log %h(aL) - L^IoiaL)) \ ah(aL)K (aL) + aI {aL)Ki(aL) 

To compute the probability p(xq) we estimate two fluxes: first at the PSD 
given by 



Jpsd = K ( udS = 2irK I — Ko(ar)rdr + 2ttk f AIo(ar)rdr 

and second the flux at the lateral boundary 

J Lat = -2irDRhu'(R) = -2tt DhC. 

Thus the probability to hit one of the receptor, where the vesicle is released 
at the center is given by 

Jpsd 
P = • 

J PSD + JLat 
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General location of vesicular release. 

To estimate the flux for a general location of vesicular release, we use the 
general expression of the Laplacian to rewrite equation ([9]), 



d 2 u(r,9) ldu{r,9) 1 d 2 u(r, 



dr 2 r dr r 2 d 2 9 " K ' //v -i u '-°J r hD 

We chose the release point on the line 9q = 0, the solution can be developed 
in a cosine series 



+ 



a 2 u{r, 0)x[q,l] 



1 



-5(r-r )5(9 



— + ^2 a n cos(n6 



(13) 



n=l 



To estimate the flux, we shall compute 



u(r,9)rdr d6 = k I I r rdrdO 



o 



o J o 



By substituting the expansion (|13p in equation(|12p and integrating with 
respect to 9 we found that ao satisfies 



d 2 a (r|r ) 1 da (r\r ) 



+ -- 

r 



dr 2 ' r dr 
We look for a solution of the form 

Alio (or) 



a 2 a (r\r )x[ ,L] 



1 



Duhr 



6(r - r ). (14) 



< r < ro 

a (r) = < A 2 /o(a^) + -82-^0(0^) r < r < L 

A 3 log(r/R) L < r < R 

When ro < L, the continuity at the point xq = (ro, 9q = 0) leads to 
Ail (ar) - A 2 I (ar ) - B 2 K (ar ) = 



(15) 



while integrating equation (|14p over (ro — e, ro + e), we obtain the condition 



d dao(r\r ) 
ro _ e dr 



a 2 ao(r|ro)rdr 



1 



r -e 



Dirh 



Taking the limit e — > 0, we get that 

A 2 ii(ar ) - -B 2 i£Ti(aro) - A 1 I 1 (ar ) 



1 



aDirhro 



(16) 
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The final condition comes from the interface r = L, where we require the 
continuity of ao and its derivative with respect to r. We get 

A 2 I Q (aL) + B 2 K {aL) - A 3 log(L/R) = 

A 2 h(aL) - B 2 K x {aL) - A 3 \ = 

aL 

From equations (|15II17|) . we obtain four independent equation for the coeffi- 
cients A\ , A 2 , B 2 , A% and the net flux can be express as 



Jpsd(xq) 



u(r\x )dS 



= k (Io(ar), 0,0,0) 



/ I (ar) -I {ar ) -K (ar ) 

-Ii(ar ) h{ar Q ) -Ki(ar ) 

I (aL) K {aL) 

h{aL) -K x {aV) 



\og{L/R) 
1 

aL 





V J 



(0,J {or), K o (ar),0) 



( I (ar) -I (ar ) -K (ar ) 

-h(ar a ) h{ar a ) —Ki(aro) 
h{aL) K a (aL) 

\ h(aL) -Ki(aL) 



\ 



In addition, the lateral flux is given by 

J Lat (x ) = -2irDRhu'(R) = -2irDhA 3 . 



log(i/i?) 



(18) 



aD 


V / 



Thus the total probability to be activated one receptor when the neurotrans- 
mitters are released at position x$ is given by 



Jpsd(xo) 



JpSd{Xq) + JLat(xo) 



(19) 



In practice, we solve (|15H17p numerically. 

When a vesicle is released at a position outside the PSD,( ro > L), ao 
has the form 

Ailo(w) < L < r 

a (r) = < A 2 log(r) + B 2 L < r < r 

A 3 log(r/R) r <r < R. 
Using similar considerations as previously, the continuity conditions lead to 
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the set of equations: 



( A x I Q (aL) - A 2 \og(L) - B2 
A ia h(aL) - ^ 
A 2 log(r ) + B 2 -A 3 log^ 








-1 

nDhro 



(20) 



Solving the linear system (|20p . we obtain the following expression for the 
flux: 



Jpsd{x ) 



(J (or), 0,0,0) 



/ h{aL) 
aI\{aL) 






- log(L) 

log(r ) 1 
-i 



1 






log a 








-1 

Dhr 



(21) 



1.2 A partial absorbing boundary condition at the PSD 

We present here our methodology to compute the partial absorbing constant 
K for an ensemble of N partially reflecting receptors of size a located on the 
PSD. When a glutamate molecule hits a single receptor, it can sometimes 
be activated or not. This condition at a single receptor is given by a partial 
absorbing condition 

-D— = K a p 
on 

where K a is the AMPA partially-reflecting activation barrier (n a = if there 
is no activation barrier and the receptor is activated upon a glutamate hit- 
ting while for K a = oo, the barrier would be so large that every glutamate 
molecule would only be reflected). The value of n a depends on the intrinsic 
properties of the AMPA binding site. 

To compute k the homogenized partial absorption coefficient, we consider 
that all the receptors are located in the PSD disk of radius Rpsd- The 
general partial absorbing boundary condition would be 

- D^- = k p on the PSD. (22) 
on 

Our criteria would be that the flux through the iV a individual receptor and 
the flux through the partially absorbing PSD should be equal. To compute 
the flux through the partial absorbing PSD, we solve the Mean First Passage 
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Time equation with the boundary condition (|22p instead of an absorbing 
boundary condition: 



DAu = -lonfi (23) 
du 
dn 



= on <90 r 



du 

—D— = k u on the PSD 

on 

To solve this equation, we use the standard method involving the Neumann- 
Green function ((37J, l38l)). In three dimensions, we find that mean first 
passage time to the PSD is approximated by 

/ x |A| f 1 ^ \ 

Thus the flux per particles is 

1 D 1 



J 



u{x) |fi| 



4i?PSD ~ l ~ 2ttkK 2 psd 



First let us consider the flux on N a AMPA receptors of size a in a disk of 
size R, which are fully absorbing. Then, the MFPT is 

|fi| N a a + f(a)R 

(25) 



ARpsdD N a a 
Thus when we equal relation (I25D with ()24p . we get the relation 

|fi| ( 1 ^ _ |0| N a a + f{a)R PS D 



D \4Rpsd 2itkpR psd J ARpsdD N a a 
which leads to the expression for the partial homogenization constant 

D N n a 



2nR psD f(a 



where f(a) = 1 - a, a = N a a 2 /R PSD . Now in general, (|30|, |40, |4l|), we 
obtain the following relation 

101 / 1 / D 

+ 



D \4RpsD N a a K a 2ita?N a 
and thus 

D 1 



^nRpSD + Ka27r a 2 N a 



(26) 
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1.3 Comparison with Brownian simulations 

Because our analytical analysis contains several approximations such as the 
averaging over the cleft height or the coefficients A\ , A2 , Bi , A3 are approxi- 
mated numerically, we evaluated the accuracy of our analysis by comparing 
the probability f|19|) with Brownian simulations (see figure [7|). We simu- 
lated Brownian particles in the same cylindrical domain as the one used 
for the analytical computation with an absorbing lateral boundary condi- 
tion. We put at the PSD a partial absorbing boundary with the condition 
—D^ = —up. At the particle level, we implemented the reflection rule (H^), 
in which particles hitting the PSD boundary are reflected with a probability 



p 



D 

and absorbed with the complementary probability, where D is the diffusion 
constant and At is the time step of the simulation. The scheme is standard 
when the glutamate molecule is inside the cleft, but when x(t)+\/2D Aw < 
at the PSD, then we use 



x(t + At) = < 



(x(t) + V2D Aw) w.p. 1 - PV At 
terminate trajectory otherwise. 



1.4 Estimation of the partial absorption rate n a using exper- 
imental data. 

Using a Markovian kinetic model for the initial binding step of a glutamate 
to an AMPAR, we estimate here the rate constant n a and the homogenized 
coefficient k with the help equation 1261 

In a two state chain model, accounting for binding and unbinding of a 
glutamate molecule to a receptor described as 

C <± O, (27) 

the forward binding rate k± is given in units of Molar In one hand, the 
binding rate is calculated by the flux formula as 

J Markov = faA^NgV' 1 f p(x)dx « kiA~ 1 N g V~ 1 TTa 2 p(x), (28) 

JdVta 



Multi-conductance AMPA receptor model 



19 



where A is the Avogadro number, p(x) is the density of glutamate near the 
receptor and V = nRh 2 is the volume of the synaptic cleft. On the other 
hand, using the diffusion model, the flux term is given 



Jdiff = N g K a / p(x)dx « NgK a Tra p(x), (29) 

JdQa 

where a is the radius of a receptor and dfl a represents the receptor surface. 
By equating the two fluxes Jdiff and J Markov, we obtain an expression for 
the partial reflecting constant 

Ka = kiA- l V~ l . 

Now, using the published value h\ = 10 7 M _1 s _1 (taken from (0)) we obtain 
that K a ~ 1.06. This two-model chain state is a good enough approximation 
even in the case where there are more states in the Markov chain. Indeed, 
the transition rates to desensitization states are lower than the open state 
so that in the short time scale, after binding (the time of interest in our 
model), the probable state is the open state. 

1.5 Mathematical details for the computation of the partial 
rate k 

We provide here the mathematical detail used in appendix for the expression 
of the partial absorbing constant k. We solve equation (|23|) asymptotically 
using the Green function: 

A x N(x,y) = -5(x-y), for x,ye9, 

dN ^y) = _J_ for XEdn ven 

du x 6 ' V G ' 

If x or y (or both) are in <9f2, then only a half of any sufficiently small ball 
about a boundary point is contained in f2, which means that the singularity 

of Neumann's function is — ; r . The Neumann's function for y £ <9f2 

2ir\x — y\ 

is given by 

where v(x, y) is a regular function. Using Green's identity and the boundary 
conditions (]22p . we obtain 

«(y)-4 / N(x,y)dx= [ N(x,y)^^-dS x +C, 
u Jn Jan 
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where 



The conservation of the flux leads (see (|39l)) to 



c= w\L u{x)dx - 



?^dS x = -B, (30) 
an a ov x D 



Following the argument in (J39|), the function N(x, y) is integrable indepen- 
dent of d£l a , whose integral is uniformly bounded, whereas C — > oo as a — > 0. 
duix^j 

Setting g(x) = — for x G dfl a and using the boundary condition (f22|) . 

dv x 

we obtain the integral equation for the flux density g(x) in d£l a , 
D 



H 



g(y) + / N(x, y)g(x) dS x = -C for y G dU a , 
Jdn a 



Using the expansion of the flux g(x) = go(x) + gi(x) + g2{x) + ■ ■ ■ , where 
gi + i(x) <C gi{x) for a — > and choose 

-2a 

go{x) = - 



cnr\ 1 



x\ 2 



a 2 



It was shown in, (|38l . 139 ) that if d£l a is a circular disk of radius a, then 



1 f 9o{x) , dSx = a forall yG ao a . 



27r f - y| 

Thus we approximate the solution by taking y at the origin, so that go(0) 
= ^ L , thus we get 

an ' ° 

The flux condition ([30]) gives 

/ g {x) dS x = -Aaa 
Jdn„ 



We conclude that 



u(x)f*C=\n\ I 7-^7 + 



27TKa 2 4Da 
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Figure Legends 

Figure Q3 

Schematic representation of the Synaptic cleft. The synaptic cleft 
geometry is approximated as a narrow cylinder of height h and the PSD is 
positioned at the center of the pre-synaptic terminal. We depicted a vesicle 
released at a distance r$ inside the Active Zone (AZ). Diffusing receptors 
can either bind an AMPA receptor or diffuse away. 
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Figure [2], 

The geometrical properties of the synaptic current. A. We decom- 
pose the synaptic current Is in a sum of 2,3 and 4 glutamate bound gluta- 
mate molecules (Is = 12+I^+Iij- hi the range of 3000-9000 glutamates, the 
contribution of each configuration is I4 > I3 > Ii- The release is from the 
center. B. The synaptic current is plotted as a function of the release dis- 
tance from the center of the synapse for one, two and three vesicles. In both 
graphs we used synapse of 500nm radius, height 30nm, PSD radius of 300nm 
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Figure [3]. 

Synaptic current (computed from[7]) versus the number of glutamate molecules. 
We show current Vs. the number of released glutamate molecules for dif- 
ferent receptor radius, a = 1.5nm, 1.8nm and 2nm . For receptor effective 
binding radius of a = 1.8nm saturation achieved at four vesicles. 
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Figure [H 

Optimal PSD radius. A. The mean current and SD are plotted as a 
function of the PSD size for three different active zones (50 nm,100 nm and 
150 nm) (the synaptic radius is 500 nm and the height is 20 nm). B, C. The 
mean number of AMPA receptors bound by 2 (resp. 4) glutamate molecules 
as a function of the PSD radius. D. Current vs. Active Zone radius. 
The PSD size is fixed at 300 nm and each curve represents 1, 2 and 3 released 
vesicles. E. CV vs. PSD size. The CV achieves a minimum when the 
PSD and the active zone are approximately equal. In that case the Active 
Zone is 100 nm and the CV minimum is achieved for a PSD of radius of 120 
nm. F. Optimal PSD radius: It is plotted as a function of the AZ radius 
obtained by minimizing the CV for a fixed AZ. 
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Figure [5J 

Synaptic current for different synapse radius.. For a fixed active zone 
radius (50 nm), we plot the synaptic current as a function of the PSD radius 
for four different sizes of synapses 200,300,400 and 500 nm (the height is 30 
nm). Doubling the size of the synapse leads to a current amplitude that 
increases from 125 to around 190 pA for a small PSD radius. (B) The 
CV is plotted as a function of the PSD radius: the CV minimum does not 
depend much on the synaptic radius. (C) The CV is plotted as a function of 
the PSD radius, for different number of released glutamate molecules. There 
is an optimal PSD size, and the CV decreases as a function of the number 
of released glutamate molecules. (D) We compare the CV as a function of 
the PSD radius curve when the number of released glutamate molecules is 
distributed according to a Gaussian distribution with a mean 3000 and a 
standard deviation of 500 (bold line) with a fixed number of 3000 glutamate 
molecules (dashed line). Fluctuation in the number of glutamate molecule 
in the vesicles does not affect much the general behavior of the CV with 
respect to PSD size. 
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Figure [6], 

A. MC Simulation For a fixed active zone radius (300 nm) a Monte-Carlo 
Brownian simulation of 500 molecules uniformly released from the AZ. Each 
experiment has 100 repetitions at different PSD size. B. CV analysis for a 
simple model. The model describes a random variable with three possible 
outcome values I\ , I2 , I3 • The probability function is given in t$5§ . The 
model describes a single receptor with three conductivity levels, depends on 
the number of bound molecules, each can bind the receptor with probability 
q. Each curve represents the CV Vs. probability of binding q with different 
values of ^3, where I\ = 1, 12 = 2 are fixed. For large values of I3, an 
optimal point emerges. The existence of optimal CV as function of binding 
probability is thus strongly connected to the nonlinear cooperative effect of 
multi-binding. 
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Figure 

Probability p(xq) to bind a receptor before exit the cleft. We com- 
pare here the analytic solution of (|19p (continuous line) with Brownian sim- 
ulations (dots) in the cylindrical domain of radius R = 0.5, height 0.02. The 
radius of the PSD is L = 0.3. 
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